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Abstract. We investigate simple random walks on graphs generated by repeated barycentric sub- 
divisions of a triangle. We use these random walks to study the diffusion on the self-similar fractal, 
the Strichartz hexacarpet, which is generated as the limit space of these graphs. We make this con- 
nection rigorous by establishing a graph isomorphism between the hexacarpet approximations and 
graphs produced by repeated Barycentric subdivisions of the triangle. This includes a discussion of 
various numerical calculations performed on the these graphs, and their implications to the diffusion 
on the limiting space. In particular, we prove that equilateral barycentric subdivisions — a metric 
space generated by replacing the metric on each 2-simplex of the subdivided triangle with that of a 
scaled Euclidean equilateral triangle — converge to a self-similar geodesic metric space of dimen- 
sion log(6)/log(2), or about 2.58. Our numerical experiments give evidence to a conjecture that 
the simple random walks on the equilateral barycentric subdivisions converge to a continuous diffu- 
sion process on the Strichartz hexacarpet corresponding to a different spectral dimension (estimated 
numerically to be about 1.74). 



The goal of this paper is to investigate the relation between simple random walks on repeated 
barycentric subdivisions of a triangle and the self-similar fractal Strichartz hexacarpet. We explore 
a graph approximation to the hexacarpet in order to establish a graph isomorphism between the 
hexacarpet and Barycentric subdivisions of the triangle. After that we discuss various numerical 
calculations performed on the approximating graphs. We prove that the equilateral barycentric sub- 
divisions converge to a self-similar geodesic metric space of dimension log(6) /log(2) ~ 2.58 but, 
at the same time, our mathematical experiments give evidence to a conjecture that the simple ran- 
dom walks converge to a continuous diffusion process on the Strichartz hexacarpet corresponding 
to the spectral dimension ps 1.74. 
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1 . Introduction and main conjectures 
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In Section [2] we develop the framework and basic results pertaining to barycentric subdivision. 
This is a standard object, intrinsic to the study of simplicial complexes, see lfT6ll (and such classics 
as Il23ll25ll28l0 . We define a metric on the set of edges of nth-iterated barycentric subdivision of 
a 2-simplex, and use this to define a new limiting self-similar metric on the standard Euclidean 
equilateral triangle. 

In Section [3} we turn to the theory of self-similar structures, as developed in |[T9l|. Using this 
theory we introduce a fractal structure, which we call the Strichartz hexacarpet, or hexacarpet for 
short. The hexacarpet is not isometrically embeddable into two dimensional Euclidean space, but 
otherwise resembles other self-similar infinitely ramified fractals with Cantor-set boundaries, such 
as the octacarpet (which is sometimes referred to as the octagasket, see (6]| and references therein), 
the Laakso spaces (see [|27ll29l and references therein), and the standard and generalized Sierpinski 
carpets (see and references therein). 

We draw a connection between the hexacarpet and barycentric subdivisions in Section|4} in which 
we prove that the approximating graphs to the hexacarpet are isomorphic to the graphs created by 
barycentric subdivisions (where the 2-simplexes of the n-times subdivided triangle become graph 
vertices, connected by a graph edge of the simplexes that share a common face). 

Section [5] discusses properties of the approximating graphs of the hexacarpet to contrast and 
illuminate connections between the hexacarpet and the limiting structure on the triangle defined in 
Section [2} In particular, we examine the growth properties of the graph distance metric. We prove 
a proposition which heuristically places the diameter (in the sense of the usual graph distance) of 
the nth level graph as somewhere between 0(2 n ) and 0(n2 n ). Our numerical analysis supports a 
conjecture for the formulas of the diameter and radius of these graphs. 

Finally in Section [6] we briefly describe numerical analysis of the spectral properties of the ap- 
proximating graphs to the hexacarpet. Primarily we calculate eigenvalues of the Laplacian matrix 
for the first 8 levels of the approximating graphs. This allows us to approximate the resistance 
scaling factor of the hexacarpet, and suggests that there is a limit resistance. We also plot approxi- 
mations to the hexacarpet in 2- and 3-dimensional eigenfunction coordinates. 

The experimental results in Section [6] strongly suggest that the simple random walks on the 
barycentric subdivisions converge to a diffusion process on K (most efficiently this can be shown 
by analyzing harmonic functions and eigenvalues and eigenf unctions of the Laplacian). Thus, our 
theoretical and numerical results support the following conjecture, which is explanationed in Sec- 
tion |6] 

Conjecture 1.1. We conjecture that 

(1) on the Strichartz hexacarpet there exists a unique self-similar local regular conservative 
Dirichlet form £ with resistance scaling factor p k-1.304 and the Laplacian scaling factor 
t = 6p; this form is a resistance form in the sense of Ki garni. 

(2) the simple random walks on the repeated barycentric subdivisions of a triangle, with the time 
renormalized by r n , converge to the diffusion process, which is the continuous symmetric 
strong Markov process corresponding to the Dirichlet form £; 

(3) this diffusion process satisfies the sub-Gaussian heat kernel estimates and elliptic and par- 
abolic Harnack inequalities, possibly with logarithmic corrections, corresponding to the 

log(6) log(Q) 
Hausdorff dimension - — — ^2.58 and the spectral dimension 2- — — rj 1.74; 

log(2) log{T) 

(4) the spectrum of the Laplacian has spectral gaps in the sense of Strichartz; 

(5) the spectral zeta function has a meromorphic continuation to C. 
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Figure 2.1. Barycentric subdivision 



We note that our data on spectral dimension is not inconsistent with the (random) geometry of the 
numerical approximations used in the theory of quantum gravity, according to the work of Ambj0rn, 
Jurkiewicz, and Loll (see LI]) at small time asymptotics: ds = 1.8 ± 0.25. This reference as well as 
[|26l use triangulations similar to those in our study to approximate quantum gravity. Therefore one 
can conclude that, with the present state of numerical experiments, fractal carpets may represent a 
plausible (although simplified) model of sample geometries for the quantum gravity. 

Acknowledgements. We are grateful to Michael Hinz and an anonymous referee for many helpful 
comments. 

2. Equilateral barycentric subdivisions and their limit 

The process of subdividing the 2-simplex using its barycentric coordinates is useful in order to 
establish an isomorphism of graphs in later sections. Information on barycentric subdivision of 
more general n-simplexes can be found in [[T6|. We adapt Hatcher's notation slightly, which is 
outlined in the following definitions. 

Consider any 2-simplex (triangle) T in the plane, defined by the vertices [vq, v 1, v 2] which do 
not all lie on a common line. The sides of T are the 1-simplexes: [v , v±], [vo, V2], [v±, V2]. 

Definition 2.1. We perform barycentric subdivision (BCS) on T as follows: First, we append the 
vertex set with the barycenters of the 1-simplexes [t> , v%), [vq, 1*2], [^1,^2] and label them 6 i, &02, &12, 
respectively. Also append the barycenter of T which is the point in the plane given by | (vq +V\ + 
v 2 ), which we denote b. Thus, is the midpoint of the segment Vj), Any 2-simplex in the 
collection of 2-simplexes formed by the set iV = {vq, v±, v 2, 601, &02, &12, b} is said to be minimal if 
its edges contain no points in iV other than its three vertices. Let B(T ) denote this collection of 
minimal 2-simplexes. Note that these six triangles are of the form bij, b] where i ^ j E {1,2,3} . 

We define the process of performing repeated barycentric subdivision on T as follows: For a 
collection C of 2-simplexes, we define B(C) = [J c€C B(c) to be the collection of minimal 2- 
simplexes obtained by performing BCS on each element of C . In this way we define the n th level 
barycentric subdivision of T inductively by B n (T ) = 5(_B n_1 (T )). 

Definition 2.2. We call the elements of B n (T ) the level n offspring of T where T is the level n 
ancestor of its 6™ offspring in B n (T ). Similarly, for any triangle T obtained from repeated BCS of 
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T , we may consider the level n offspring ofT to be the collection B n (T). We use the terms child, 
(resp. grandchild) to denote the level 1 (resp. level 2) offspring of T. Likewise, we use the terms 
parent, (resp. grandparent) to denote the level 1 (resp. level 2) ancestor of T. We will use t C T 
to denote that t is a child of T, and when necessary t C T C V to denote that t is a child of T and 
a grandchild of T'. If s and t are both children of T, then we say that s and i are siblings. 

Definition 2.3. For any triangle T = [a, 6, c] we define the boundary of T to be the union of its 
sides, which we denote dT = [a, b] U [6, c] U [a, c] . A level k offspring t of T is said to be a boundary 
triangle for T or on the boundary of T if a side of i lies on <9T. For a given triangle T, we say that a 
level fc offspring of T is special with respect to T if it is on the boundary of T and contains a vertex 
of T. Note that all 1 st level offspring are special with respect to T, and that every special offspring 
has exactly on special offspring. 

Definition 2.4. We say that two level n triangles are adjacent if they share a side. Given a level 
n triangle T = [vq,v\, v 2 ], we know the children of T are of the form [vi, bij,b] where i ^ j 6 
{1, 2, 3}. We say that two children of T are vertex adjacent if their common side is a segment 
connecting the barycenter of T to one of the original vertices of T ([vi, b] for some i G {1,2,3}) 
and that two children of T are ^z'Je adjacent if their common side is a segment connecting the 
barycenter of T to one of the barycenters of the sides of T ([6^, 6] for some i G {1,2, 3}). 

Note that each application of BCS on any triangle T produces six new offspring, so we have 
| B(C) | = 6- \C\ for any collection of 2-simplexes C. Similarly, starting with T we have inductively 
that |_B n (T )| = 6 n . Thus, there are 6 n level n offspring of T . Also, we see that each triangle in 
B n (T ) is adjacent to at most three other triangles in B n (T ) and that if t G B n (T) is not on the 
boundary of T, then t is adjacent to exactly three other members of B n (T). On the other hand, if 
t G B n (T) is on the boundary of T, then t is adjacent to exactly two other members of B n (T), 
namely its vertex adjacent sibling and side adjacent sibling, and possibly one other triangle not in 
B n {T] ) but adjacent to T. 

Proposition 2.5. The following are immediate consequences of the above definitions. 

(1) The number of level n boundary triangles ofT is 6 • 2 ra_1 . 

(2) The number of adjacencies among B n (T ) is §[2 • (6 • 2"" 1 ) + 3 • (6" - 6 • 2™- 1 )] = 
2 „_i _ ^ n+ i _ 

(3) If s C S and t C T are adjacent, then either S = T or S is adjacent to T. 

(4) Ift C T, then exactly one child oft is special with respect to T. 

Theorem 2.6. On the triangle T there exists a unique geodesic distance d^x, y) such that each 
edge of each triangle in the subdivision B n (T ) is a geodesic of length 2~ n . The metric space 
(T , doo(x, y)) has the following properties: 

(1) (T , doo{x, y)) is a compact metric space homeomorphic to T with the usual Euclidean 
metric \x — y\; 

(2) the distance d^x, y) from any vertex of a triangle in _B„(T ) to any point on the opposite 
side of that triangle is 2~ n ; 

(3) there are infinitely many geodesies between any two distinct points; 

(4) (T , rfoo(x, y)) is a self-similar set build with 6 contracting similitudes with contracting 
ratios \; 

(5) The Hausdorff and self-similarity dimensions of(T , d^x, y)) are equal to j^f)- 
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Proof. We choose a particular triangle T = [v , v±, v 2 ] and construct a metric doo on T by approx- 
imating with metrics d n on dB n (T ), the union of edges of the triangles in B n (T ). We have that 
dB n (T ) is a 1-simplicial complex, a quantum graph, and a one dimensional manifold with junction 
points at the vertexes of B n (T ). There is a unique geodesic metric d n (x, y) on dB n (T ) such that 
the length of each edge is 2~ n . It is easy to see that for any x,y E dB n (T ) and any k > we have 
d n (x,y) ^ d n+ k(x,y). Moreover, by induction one can show that, if x,y vertices of B n (T ) and 
k > 0, then we have the compatibility condition d n (x,y) = d n+ k(x,y). Therefore on the union 
U^ =0 dB n (T ) there is a unique geodesic metric d^x, y) such that each edge of each triangle in the 
subdivision B n (T ) is a geodesic of length 2~ n . 

Since the diameter, in the Euclidean metric, of 2-simplexes of r-times repeated barycentric sub- 
divided equilateral triangle is bounded by (2/3) r , U^ =0 dB n (T ) is dense in T with respect to the 
Euclidean topology (see Chapter 2, Section 1 of IfToTO . Thus d^ can be extended to T by continuity. 

Next we notice that, with the metric d^, T is self-similar in that if T = [a , ai, a 2 ] is a 2- 
simplex of B n (T ), then the map of simplexes / sending Vi — >■ a«, and extending by linearity, is a 
homeomorphism with (x, = 2~ n d O0 (f~ 1 (x),f~ 1 (y)). 

To prove (2), we note that it is enough to show that for, x a corner of T , and y a vertex of B n (T ) 
on the side of T opposite x, d oc (x,y) = 1. This is sufficient because of self-similarity and that 
such vertexes become dense in the side of T which is opposite x. Any such y can be connected to a 
barycenter of a 2-simplex of _B n _ 1 (T ), call this barycenter by an edge of length 2~ n . In turn, 
x„_i can be connected to x n _ 2 , a barycenter of a 2-simplex of B n _ 2 (T ) by an edge of length 2~ n . 
In turn x n _ 2 can be connected to a barycenter of a 2-simplex of _B n _ 3 (T ), called x„_ 3 , by a path 
of length 2~ n ~ 1 . Proceeding by induction we get that x can be connected to y by a path of length 

2 _1 + 2~ 2 H h 2~ n+1 + 2~ n + 2" n = 1. 

This is precisely d^ (x, y) because it can be shown that any path from x to y on dB n (T ) has to pass 
through at least 2™ — 1 vertexes. This argument implies, in particular, that the diameter of T in the 
metric d^ is equal to 1, and so the metric is finite even when continued from U^L <9-B n (T ) to T . 

The same kind of argument also proves (3). If y is a vertex of B n (T ), then it is a vertex of 
B n+k (T ) for all k ^ 0, and the collection of paths described above differ depending on our choice 
of n and k. It is then easy to construct infinitely many geodesies with the same length between 
arbitrary distinct points of T . 

To prove (1), we use the insight into the properties of d^ described above. In particular, the d^ 
ball of radius 2~ n centered at x, a vertex of B n (T ), is contained in the union of all 2-simplexes 
of B n (T ), which contain a; as a vertex. The union of these triangles clearly contains a Euclidean 
open set containing x, and is contained in the Euclidean ball around x of radius (2/3) n . Since 
U^i 1 S Tl (T ) is dense in both metrics, this proves that the metrics are equivalent in the sense that 
they induce the same topology. This proves (1), and has the added bonus of proving that T is 
compact with respect to d^, and thus d^ is complete. 

The compactness of T with respect to d^ proves (4) because T is the union of the 6 contraction 
maps from T onto each 2-simplex of Bi(T ), and so it is a self-similar structure (see the next 
section). Each of these maps has contraction ratio 1/2, and thus using a calculation found in |flQ|. 
one discovers that the self- similarity and Hausdorff dimension of T is log(6) / log(2), which proves 
(5). " ' ' □ 
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3. Self- Similar structures and the Stricharz hexacarpet 

First we introduce some notation that we will use. We denote X = {0, 1, . . . , N — 1}, called an 
alphabet, and 

X n = {X1X2 ■ ■ ■ x n I Xi G X} 

will be the words of length n. Also, we take 

oo oo 

X* = |JX n and, E = Jjx 

n=0 i=l 

Naturally, the set X has discrete topology, and E is given the product topology (i.e. the topology 
whose basis is sets of the form YYiLi A-u such that A{ = X for i > M for some M). There is even a 
natural metric on E, defined below. 

Proposition. Fix a number r G (0, 1). For w = W1W2 ■ ■ ■ and v = V\V2 ■ ■ ■ in E, we define 
S r (w, v) = r n where n = min {£ : W£ ^ Vg} with the convention that 8 r (w, w) = 0. Then 5 r is a 
metric on E. Additionally, the maps <Ji{w) = iw, for i G X, is a contraction with Lipschitz constant 
r. 

This is proven as Theorem 1.2.2 in [[191. The work [fT9l introduces the the theory of self-similar 
structures and is developed in the context of contractions on metric spaces. We also use the defi- 
nition of self- similar structure set forth in the above paper. In the rest of this section, we shall take 
S = 5i/2, to be the metric that makes o,i contractions with Lipschitz constants 1/2. 

Proposition 3.1. Let K be a compact metrizable space, let X be a finite indexing set for Fj : K — >• 
K continuous injections such that K = U i( zxFi(K). We call the triplet C = (K,X, {-Fj} igX ) a 
self-similar structure on K if there is a continuous surjection n : E — >■ K such that the relation 
Fj o 7T = 7T o <7i, where di{w) = iw, holds for all w G X*. 

We define the n th level cells of K, K w = F W (K) for w = wiw 2 ■ ■ ■ w n G X n , where F w = 
F Wl o F W2 o • • • o F Wn . In particular, if K is some quotient space of E, where n is the quotient map 
where (Tj's are constant on the fibres of7rforalH G X, then we can define Fi = n o OiOix _1 , creating 
a self-similar structure. It is in this way we shall define a fractal. 

We define the equivalence relation ~ by the following relations: Let X = {0,1,...,5} where x 
is any element in X* and v G {0, 5} w . Suppose i,j G X and j = i + 1 mod 6. Then, if i is odd, 

(3.1) xi3v ~ xj3v and xiAv ~ xj4v. 
If i is even (j is still i + 1 mod 6), then 

(3.2) xilv ~ xjlv and xi2v ~ xj2v. 
We define K := E/ ~. 

We may also define K in an alternate way, which we shall call K. The equivalence relation on 
K is defined by 

xiy ~ xjz 

for x G X*, i, j G X, and z, y G E, where 

j — i + 1 mod 6 
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Figure 3.1. G 2 generated by ( |3.1| ), p.2| ) on left, and ( |3.3| ), ( |3.4[ ) on right 
and, if i is odd, then 

(3.3) yk = i + 2ori + 3 mod 6 and 
and if i is even, then 

(3.4) yk = i + lori + 2 mod 6 and 



1 mod 6 if Dk = i + 2 mod 6 

2 mod 6 if = 2 + 3 mod 6. 



2/c 



if t/fe = « + l mod 6 
1 mod 6 if = i + 2 mod 6. 



Proposition 3.2. 77?e space (K, X, {o"j}) z's a self-similar structure, where ai : K ^ K is defined, if 
7T : £ — >■ if w projection associated with the equivalence relation ~, o"i(7r(x)) = 7r(<jj(a;)). 

Proof. Since x ^ y implies that crj(x) ~ cr^y) for all x, ?/ G S and i E X (this can be seen in the 
definition of ~), cr; : K — > K is well defined. The only thing left to check is that K is metrizable. 
This follows because we can define a metric 5 on K by 

S([x],[y])= M {5 r (x,y)}. 
This metric is well defined because the equivalence classes of ~ contain at most 2 elements. □ 



Proposition 3.3. The equivalences defined in equations p.l[ ),( [3~!2] ) and equation ( |3.3[ ), ( |3.4| ) provide 
two definitions of K which are equivalent. 

Proof. We show that there is a self homeomorphism /:£—>•£ which transforms the equivalence 
in 



3.1 and 3.2 into the equivalence 3.3 and 3.4 This map / is given by f(u) = v, where 



vi := ui and v m := (-l) am - x u m + a m _i, 



where a fe = ^7=1 u j- This map is continuous, since how / acts on the nth letter is independent of 
any future letters, 6(x, y) < (f(x), f(y))- In fact, since / acts bijectively on X n (which is easy to 
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check by induction), this means that / is an isometry with respect to 5 r , and is thus a bijection as a 
map of from £ to itself. Moreover, / also descends to an isometry with respect to 5, but to see this 
we need to know that / preserves the equivalence relation, which we will show shortly. 

Another way of seeing that / is bijective comes from that fact that we can define an inverse 
f~ x (v) = u where 

Mi := v t and u m := (-l)° !m - 1 (t> m - a m _i). 
Showing that the two equivalences produce identical quotient spaces is a matter of showing that / 
and f^ 1 preserve equivalence, i.e. u ~ v with respect to 3.1 3.2 implies that f(u) ~ f(v) with 
respect to 3.3|3.4 This involves checking a few cases, we provide an example of such a case. 

For the remainder of the proof, all equalities are assumed to be mod 6. Now suppose that vP^ 
is of the form xi3w and is of the form xj3w where the j — % + 1 and i is odd, as in 3.1 and 
we assume that the i and j are in the nth position. Finally we assume, letting = Y^i=i anc ^ 



a 



iP = Y^=i mat a n~i = a n-i i s °dd. We need to vary all of the above for various cases of 
equivalence. 

Let = f(u®) for I = 1, 2, then for k < n, = vf \ and vh = Vn^ — 1. Furthermore, Vn^ 
and is odd, so t^+x = + 3, and 

v { nli = + 3 = - (j - 1) + 3 = + 4 = - 2 mod 6, 

this is consistent with 13 .31 

For k > n + 1, if we have = uf^ = 0, then v% = u^ x and v k = vjpi- The first instance 

(if it happens at all) where ujj. = u% = 5, aj^ = + 3 will be odd, so u£ = t4 + 2 and 
a k-i = a ^ + ^ will be even, so vj*' = — 1. Any further instance where u^' = u£' = 5, with 
add or subtract 1 to alternatively and any instance where — u ^ — will add nothing. In this 
way, we have = Vn + 2 or + 3 with v% = VrP — 1 or — 2 respectively. This shows 



that ~ v^ 2 > in the sense of 3.3 



It is in this way that all cases can be checked. □ 

We define the shift map a : E — > S by a(xiX2X3 ■ ■ ■ ) = X2X3 ■ ■ ■ . In this way, the concatenation 
maps (Tj can be seen as branches of the inverse of a. 



Lemma 3.4. If we use the definition of K from equation 3.1 and 3.2 the shift map a descends to a 
well defined map from K to K, which we also call a. 

The proof of the above is a matter of showing, for x, y E £ that x ~ y if and only if a(x) ~ cr(y). 



This is easily verified to be true in the equivalences in 3.1 and 3.2 



We turn to topological properties of the space, a natural question is what does the intersection of 
two neighboring cells look like. 

Proposition 3.5. Ifi, j e X and i 7^ j, then either <Ji{K) D <7j(K) is empty or homeomorphic to 
the middle third Cantor set. 

Proof. The set <7i(K) consists of infinite words beginning with the letter i, the intersection with 
<Jj(K) is the set of words which begin with an i which are equivalent to the words which start with 
a j. If i ^ j ± 1, then this intersection is empty, Since there is no loss in generality, we assume that 
j — i + 1, if we further assume that i is odd, then 3.1 tells us that Ui(K) fl crj(K) is given by the 
words 

i3v ~ j3v or iAv ~ jAv 



where v G £ is an infinite word consisting of O's and 5's. By ignoring the leading i or j, this set 
is naturally homeomorphic to the shift space of {0, 1} — which is in turn homeomorphic to the 
middle third Cantor set (see, for example, IPT910 . The case where i is even follows the exact same 
argument. □ 

If we examine planar realizations of the approximating graphs, we see a large "hole" in the center. 
We see that this hole consists of truncation of the elements of the set C consisting of words of the 
form ijv where i, j G X where i is any letter and j is a 3 or 2, and v G £ is a word consisting of O's 
and l's. 

Proposition 3.6. The set C is homeomorphic to the circle. 

Proof. If we consider C fl <Ji(K) (the subset starting in i), then we have that «210 ~ «310 and 
zx3010 ~ ix3llQ. for any finite word x. The shift space of {0, 1} with xOlO ~ xllO is homeo- 
morphic to the unit interval (this is seen in the limit space of the Grigorchuk group, see [|24| section 
3.5.3). 

In this homeomorphism, the endpoints of the interval are and 10. So we have two copies of 
the unit interval corresponding to the sets starting with i2 and i3, which are identified at one of the 
endpoints. This shows that each set C fl <Ji(K) is itself isomorphic to the interval. These intervals 
are in turn identified at their endpoints, z20 ~ k20 and i30 ~ £30, where k = i ± 1 and £ = i^f 1, 
depending on whether i is odd or even. □ 

From the results above we obtain, in particular, the following theorem. 

Theorem 3.7. The self-similar set K, defined above and called the Stichartz hexacarpet, is an 
infinitely ramified fractal not homeomorphic to T . 



Proof. This follows by proposition 3.5 Since K can be disconnected into arbitrarily small pieces 
by removing a topological Cantor set, it must be topologically one dimensional. On the other hand, 
T is topologically 2 dimensional, thus the spaces cannot be homeomorphic. The general theory of 
topological dimension can be found in [fTTTl (in particular see Definitions 1.1.1 and 1.6.7, Theorems 
1.4.5, 4.1.4 and 4.1.5). □ 

4. Graph Approximations and Isomorphism 

In this section we show that the self-similar structure from the last section can be approximated 
by graphs constructed from repeated barycentric subdivision of a 2-simplex. 

We make this precise by constructing approximating graphs to the hexacarpet K. Taking X = 
{0,1,..., 5} as in the last section, we define the graph G n = (X n , E n ) with X" as the set of vertices. 
We define the edge relations E n where two words u and w are connected if there are x, y G £ such 
that the concatenated words wx ~ uy according to the equivalence defining K, as in equations 
p.l[ ) and p.2[ ). We can alternatively define the the set of vertices to be the set of n th level cells 
of K, where (K w , K u ) are in the edge relation if K w fl K u ^ 0. In this way, we can think of the 
vertices of G n as being nth level cells of K. 

We now exhibit partitions of the vertex set X n and edge set E n which will be useful in discussing 
edge relations: Let Wx = {x\X 2 ■ ■ ■ x n \ Xi G {0, 5} , 2 < i < n} be the set of words of length n 
whose second through last letters are O's or 5's. For each x = X\X 2 • • • x n G X n \W 1 there is at 
least one i G {2, 3, • • • , n} such that {0, 5}. So we may define the function / : X"\Wi — >• 
{2, 3, • • • , n} by I : x±x 2 ■ ■ ■ x n — >■ max {2 < i < n \ Xi ^ {0, 5}}. Now for 2 < k < n define 
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Wfc = / _1 (A;) = {X1X2 ■ ■ ■ x n | Xk ^ {0, 5} , Xi E {0, 5} , k < i < n}. These are the words which 
end in exactly (n — k) 0s and 5s. 

From the equivalence relation ~ defined on £ we recover the following edge relations on X" by 
truncating the relations after the n th coordinate: 

Fi = {{xi, xj} | x E X n_1 , j — i + 1 mod 6}. 

Ffc = , xjav } \ x E X k ~ 2 , j — i+1 mod 6, a = 3 or 4, v E {0, 5} n ~ fe j, for 2 < 

k < n. 

We simplify the edge relations by writing {xiav, xjav}. Here a E {3, 4} is the same in both 
components of the relation when i odd, j = i + 1 mod 6. We take a E {1, 2} to be the same on 
both sides when i even, j = i + 1 mod 6. 

We now collect some information about the cardinalities of the vertex and edge sets. 

Proposition 4.1. The following are apparent from our construction: 

(1) |Wi| = 3 ■ 2 n and \\N k \ = Q k ~ l ■ 4 ■ 2 n ~ k = 3*- 1 ■ 2 n+1 ,for 2 < k < n. 

(2) |Fi| =6 n and\F k \ = 6 fe-1 ■ 2 n ~ k+1 ,for 2 < k < n. 

(3) |E n | =ELiN = ■ (3 m+1 - 3). 

Proposition 4.2. Every vertex in Wi /zas degree two, with both edge relations in Fi. For k > 2, 

every vertex in \N k has degree three, with two edge relations in Fi anJ ?/z?>J m F^. 

Proof. First we note that each vertex has exactly two edge relations in F 1? namely 

{xix 2 ■ ■ ■ x n -iy, x} and {x, X\X 2 ■ ■ ■ x n ^\z) where y = x n — 1 mod 6 and z = x n + 1 mod 6. 
In addition, for all 2 < k < n each vertex in W^. has one additional edge relation in F k - By 
construction, the vertices in Wi do not have any edge relations in F k for all k ^ 1. □ 

We now collect our information on the graph approximation of the hexacarpet and the graph 
constructed using repeated Barycentric subdivision of the 2-simplex in order to establish an iso- 
morphism between them. We begin by introducing some information which will be useful in the 



proof of Theorem 4.7 at the end of this section. 



Proposition 4.3. There exists a labeling ofB n (T ) with the strings in X™ that establishes a bijection 
between the two sets. 

Proof. Label B(T ) with the elements in the alphabet X = {0, 1, 2, 3,4, 5} cyclically so that we 
have the edge adjacencies {1,2}, {3, 4} , {5, 0} and the vertex adjacencies {0, 1} , {2, 3} , {4, 5}. 
We call this construction a standard labeling of the offspring of T with the letters of the alphabet X. 
We note that there are six standard labelings of the offspring of T , and each is uniquely determined 
by labeling any one child. 

We choose an arbitrary child of T to be labeled and construct a standard labeling of the re- 
maining children. Thus we have labeled the triangles of B(T ) bijectively with the words in X 1 via 
the standard labeling map $i : B(T ) — > X. For n > 2, we shall define an inductive labeling of 
B n (T ) with the words in X n and establish the bijection $ n : ^"(To) -)■ X n . 

By assumption, for each triangle t E B"^ 1 (T ) we have an associated unique word, x = X\X 2 • • • x n _\ 
inX™" 1 ^. $ n _i : B n ^ 1 (To) — > X n_1 is a bijection). We will label the offspring of t with the words 
xO, xl,x2, x3, xA, x5 as follows (where xO denotes the word of length n, X\X 2 • • • x„_i0 E X n ): 

Let T be the parent of t. From above we know that exactly one child of t is special with respect 
to T. We assign this triangle the word xQ and label the other children of t according to the standard 
labeling fixed by xO. Therefore, to each element of B n (T ) we have associated a word of X n . 
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To show that this is an injective labeling, assume that there are two level n offspring s and t of T 
which have the same label, say x\x 2 ■ ■ ■ x n G X n . By the induction assumption, there is exactly one 
triangle T G -B n_1 (T ) with the labeling x±x 2 ■ ■ ■ £ n -i- This means that s and t are both children 
of T which have the same label x n by the standard labeling of the children of T. Thus, s and t are 
the same triangle. So we have $ n is an injective map between the finite sets B n (T ) and X n . Since 
\B n (T )\ = \X n \ = Q n , we see that $ n is a bijection as desired. □ 

Definition 4.4. For any triangle s G B n (T ) there exists a unique chain of ancestors, s = s n C 
s n -i C • • • C Si C So = T called the family tree of s. 

Proposition 4.5. If s, t G B n (T) and s is adjacent to t then there is some maximal < m < n 
such that s m = t m in the family trees for s and t (i.e. s and t are both level £ offspring of s m = t m .) 
In particular, as k level triangles, s k and t k are equal for all < k < m and are adjacent for all 

m < k < n. 

Proof. For any adjacent s and t we have that s = t = T so this assignment is well defined. Since 
each triangle has a unique parent, if s m = t m for some < m < n, then s m _i = t m -i. Thus, we 
have s k and t k are equal for all < k < m. To see that s k and t k are adjacent for all m < k < n we 
note that if s k C s k _i and t k C t k _i are adjacent, then s fc _x and t k -i are either adjacent or equal. 
For all n > k > m + 2 we have s fc -i and are not equal by assumption. Also by assumption, 
s n and t n are adjacent. Therefore, by induction we see that the final statement of the proposition 
holds. □ 

Lemma 4.6. Let s and t be adjacent level n offspring ofT with s C S C S' and t C T C T'. 
Then exactly one of the following is true of the labeling of s and t. 

(1) If S = T has the label x e X ra_1 , then, either s has label xi and t has label xj for some 
ieX,j — i + l mod 6, or t has label xi and s has label xj. So the addresses of s and t 
differ only in their last letter by 1 mod 6. 

(2) If S 7^ T and S' = T' has the labeling x e X n ~ 2 , then s is labelled xia and t is labelled 
xja. Ifi is even, then a is either a 1 or 2 (in the addresses of both s andt) and ifi is odd, 
then a is either a 3 or 4 (in both addresses). 

(3) IfSj^T and S' ^ T', then the addresses for s and t end in the same letter, which is either 
aO or 5. 

Proof. The proof of (1) is immediate from the standard labeling procedure. We have forced the 
children of every triangle to observe the adjacencies {i, j} for all j — % + 1 mod 6. 

For (2), we see that S and T satisfy the hypothesis of (1), therefore without loss of generality we 
may assume S = xi and T = xj for some j — i + 1 mod 6. Now i is either even or odd. If i is 
even, then we have S and T are vertex adjacent siblings (this was also forced by our construction). 
Thus, their common edge is of the form u = [v k , b], k e {0, 1, 2}, where v k (resp. b) is a vertex 
(resp. the barycenter) of S' = V . By the standard labeling of the offspring of S and T, we see that 
the labels of s and t end in either a 1 or 2. Assume for contradiction that the last coordinate of s 
and t are different. Thus, without loss of generality, we have that s = xil and t = xj2. We see that 
s = [v k , bs, b u ] and t = [b, 6 T , b u ] where b u is the barycenter of the segment u and b$ (resp. br) is 
the barycenter of S (resp. T). Since s and t have only one vertex in common, it is impossible for 
s and t to be adjacent. Now we see that the address for s is xia and the address for t is xja where 
a G {1, 2} is the same in both addresses. If % is odd, then S and T are side adjacent siblings and 
the proof follows as in the case when % is even. 

n 



Figure 4.1. The labeling B 2 (T ) (left) and the graph isomorphism to G 2 



For (3), we know that 5" and V are adjacent, so without loss of generality, let their common side 
be [vq, Vi). After BCS we see that S and T must have either [t> , 601] or [v\, 601] as their common side. 
Again without loss of generality, assume that the common side of S and T is [vq, 601] . Subdividing 
further we see that the common side of s and t must be either [t> , m) or [m, 6 i]> where m is the 
barycenter of [v , 6 i]- Since [t> , m) C [v , v\], if [v , m] is a side of s and t we see that s and t 
are each special with respect to their respective grandparent, S' and T". Thus, the last letter in the 
addresses of both s and t must be by construction. On the other hand, if [m, b i] is the common 
side of s and t, we see that s and t are both side adjacent to their siblings who received a label of 0. 
Thus, the last letter in the addresses of s and t must be 5 by the standard labeling construction. □ 

We are now in a position to define the desired isomorphism of graphs. 

Theorem 4.7. Let B n (T ) be the vertex set of a graph where two level n offspring ofT are con- 
nected by an edge if and only if they are adjacent as level n offspring ofT . This graph is isomorphic 



to G n with the isomorphism given by Proposition 4.3 



Proof. We already have that $„ : B n (T ) — > X n is a bijection between the vertex sets of the two 
graphs. It remains to show that $„ preserves the adjacency structure of the two graphs. Most of 
the hard work was done in Lemma 4.6[ We must now show in particular that if s, t G B n (T ) are 



adjacent triangles, then $ n (s) and $ n (t) satisfy an edge relation in F k for some 1 < k < n. 

Let s = s n C s n -i C • • • C si C s = T and t = t n C t„_i C • • • C ti C t = T be 
the family trees for s and t, respectively, and let m be maximal with respect to s m = t m . Assume 
^m(sm) = & m (tm) = Xi%2 ' ' ' % m £ X m . First suppose that m = n — 1 and $ n _i(s n _i) = 



^n-i(^n-i) — X1X2 ■ ■ ■ x n -i G X™ . Then, by part (1) of Lemma 4.6 without loss of generality 
we have $ n (s) = X1X2 ■ ■ ■ x n -\i and $„(t) = x\x<i ■ ■ ■ x n -ij where j = i + 1 mod 6. Thus, 
{$ n (s), $„(*))} G F l5 as desired. 
Now suppose that m = n — 2 and $ n _2( s n-2) = & n -2(s n -2) = x E X" -2 . We apply part 



(1) of Lemma 4.6 to s n _i and t n _i to obtain $„-i(s n _ 1 ) = and $„-i(t n _i) = xj, where 



j = i + 1 mod 6 and i is either even or odd. We apply part (2) of Lemma |4T6| to s n and t n to see 
that $ n (s) = ■ • • x n -2ict and $ n (t) = X1X2 • ■ ■ x n _2ja, where j — i + 1 mod 6, a is the 
same in the addresses of both s and t, and a as above. Thus, {$„ (s), $ n (t))} e F n , as desired. 
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Finally, suppose that m < n—3. Fromparts(l)and(2)ofLemma4.6 we know that $ m+2 (s m +2) = 
xia and $ m +2(^m+2) = xja, where $ m (s m ) = $ m (t m ) = x G X m , j = i + 1 mod 6, and a 



is as above. For 3 < k < n 



m. 



we see that s m+ t and t m+ k satisfy the conditions of part (3) of 



Lemma 4.6 



Thus, the last label in the addresses of both s m+ k and t m+k is either or 5. Induc- 
tively we have $„(s) = xiav and $„(£) = xjav , where j = i + 1 mod 6, a is as above, and 
v G {0,5} n " m " 2 . Thus, <£>„(£))} G F m+2 , as desired. 

We now see that to each edge in £> n (T ) there corresponds an edge in E n . We verify from Propo- 
sitions |2.5| and |4.1| that the number of edges in each graph is the same. Therefore, we have an 
isomorphism of graphs given by $ n . □ 

5. Graph Distances 



Each G n inherits a natural planar embedding from B n (T ), see figure 4.1 Two interesting fea 



tures of these embeddings are their central hole and outer border. Figure 5.1 generated by the 
computer program Mathematica, shows a deformed embedding, accentuating these features. In 



this section we use the isomorphism established in Theorem 4.7 to derive formulas for the length 
of paths that follow the outer border and the central hole. We call these paths outer and inner 
circumference paths. 




Figure 5.1. Graph G 4 , with typical radius (left) and diameter (right) paths highlighted. 



Definition 5.1. Let Y n denote the collection of level n boundary triangles of T defined in Defini- 
tion 2.1 Define the extended level n boundary triangles ofT as the set of all t G B n (T ) such 
that t intersects dT in exactly one point. We denote the collection of extended level n boundary 
triangles of T by Z n . Define the outer circumference path of T , Out n , to be the cycle that crosses 
each triangle in Y n and Z n . Define the inner circumference path of T , Inn n , to be the cycle that 
crosses each triangle containing the barycenter b of T as a vertex. 

Proposition 5.2. We have the following formulas for the length of Out n and Inn n : 
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(1) \Out n \ = 3n-2 n 

(2) \Inn n \ =3-2" 

Proof. We count inductively the length of the outer circumference of T and inner circumference 
of T . The statements about G„ follow from the isomorphism. When we subdivide any triangle, 
exactly two of its children contain any one vertex of the parent. In particular, we have that for each 
level n triangle in Z n exactly two of its children are in Z n+ i. Similarly, for each level n triangle 
containing b exactly two of its children contain b. Since \Inn\\ = 6, it follows inductively that 
\Inn n \ = 6 • 2"' 1 = 3 • 2 n . 
Note that when we subdivide any triangle in Y n , exactly two of its children are in Y n+ i and exactly 



n 



two more of its children are in Z n+ \. From Proposition 2.5 we know that \Y n \ = 6 • 2 n =3-2 
thus 

n-1 

\Z n \ =2|y n _ 1 | + 2|Z„_ 1 | =£>'|^_!| = (n- 1)3-2". 

So \Out n \ = \Y n \ + \Z n \ = 3n-2 n . □ 

Next we recall some standard definitions regarding distance on an arbitrary finite graph G = 
(V, E), see [USEE] for references. 

Definition 5.3. The graph or geodesic distance d(x, y) between two vertices x, y G V is the length 
of the shortest edge path connecting them. The eccentricity £(x) of a vertex x G V is defined as 
£(x) := max{d(x, y) : y G V}. Now the diameter and radius of a finite graph G can be defined 
as V(G) := m&x{£(x) : x G V} and TZ(G) := min{£(x) : x G V}, respectively. A vertex x G V 
is called central if £{x) = 71(G) and peripheral if £(x) = V(G). If the length of the shortest 
edge path connecting a central vertex to another vertex equals the radius of a graph, then we call 
that path a radius path. A diameter path is defined analogously. Note that a radius or diameter path 
connecting two vertices need not be unique. 

Using Mathematica's graph utilities package we were able to compute the radius and diameter 
of G n for 1 < n < 9 and to plot radius and diameter paths. See figure 5.1 for some examples. On 
observing these paths we noticed that the typical radius path of G n+ i is composed of a partial path 
along Inn n+ i and an approximate diameter path of G n . This suggests the equation 

(5.1) 7Z(G n+1 ) = \PInn n+1 \ + V(G n ) - Dad Jn 

where \PInn n+ \ \ is the length of the partial path along Inn n+ \ and Dadj n is the adjustment needed 
for agreement with our data. Similarly, the well known fact that 71(G) < V(G) < 271(G) for any 
finite graph G suggests the equation 

(5.2) V(G n ) = 2K(G n )-Radj n 

where Radj n is the adjustment needed for agreement with our data. Both adjustments and | PInn n+ i \ 
appear to obey simple recurrence relations, see tables 5.1 and 5.2. Solving these relations using the 
standard techniques gives the likely formulas 

(5.3) \PInn n+1 \ = 2 n+1 + 1, 

(5.4) Dadj n = l(2 n -(-l) n -3), 

o 
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(5.5) Rad Jn = 1(7 ■ 2" + 2(-l) n + 6). 

o 

Combining equations (5.1) through (5.5) yields the recurrence relation 

(5.6) n{G n+l ) = 2K{G n ) + \{2 n+2 - (-1)" + 3). 

o 

Solving this recurrence relation results in explicit formulas for lZ(G n ) and V(G n ). 

Conjecture 5.4. We conjecture the following formulas for the radius and diameter of G n : 

(1) K(G n ) = ^(2 n+1 (13 + 3n) + - 9) 

(2) V(G n ) = |(2^- 1 (31 + 12n) + 2(-l)"- 1 - 18) 

Remark 5.5. If conjecture 5.4 is true then the following are easily seen to hold: 

(1) Asn -> oo, \Inn n \ = o{V(G n )) yet \Out n \ = 0(V{G n )) 

(2 ) lim = ? 

n^oo T>(G n ) 2 

(5; ft(G n+2 ) = 4^(G n+ i) - m(G n ) - |(1 - (-1)") 
W P(G n+2 ) = 4V(G n+1 ) - 4V(G n ) - 2(1 + (-1)™) 
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42 
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= 2 n+1 + 1 


+ 
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V{G n ) 


- I(2»-(-l)»-3) 



Table 5.1. Observed relation between 1Z(G n+ i), \PInn n+1 \, and V{G n ). 



6. Numerical Data, Concluding remarks and conjectures 



From Theorem 4.7 , the level n hexacarpet G n and the nth Barycentric offspring of a triangle are 
isomorphic. The following results assume that we are working with vertices from G n , but without 
loss of generality, we can assume they are cells of B n (T ). 

We look to solve the eigenvalue on the hexacarpet on G n 

(6.1) — A n u(x) = Xu(x) 

at every vertex in G n . For a finite graph, G n , the graph Laplacian — A n u(x) is given as 

(6.2) - A n u(x) = ^2(u(x) - u{y)) 

n 

15 



Tl 
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= 4550 
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V(G n ) - 


= 2TZ(G n ) 


- \{7 ■ 2 n + 2(-l) n + 6) 



Table 5.2. Observed relation between V{G n ) and lZ(G n ). 
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71 = 7 


n = 8 


1 


0.0000 


0.0000 
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1.0000 


1.0000 


3 


1.0000 


1.0000 


4 


3.2798 


3.2798 


5 


3.2798 


3.2798 


6 


5.2033 


5.2032 


7 


7.8389 


7.8386 


8 


7.8389 


7.8386 


9 


8.9141 


8.9139 


10 


8.9141 


8.9139 


11 


9.4951 


9.4950 


12 


9.4952 


9.4950 


13 


17.5332 


17.5326 


14 


17.5332 


17.5327 


15 


17.6373 


17.6366 


16 


17.6373 


17.6366 


17 


19.8610 


19.8607 


18 


21.7893 


21.7882 


19 


25.7111 


25.7089 


20 


25.7112 


25.7091 



Table 6.1. Hexacarpet renormalized eigenvalues at levels n = 7 and n = 8. 



for every vertex that neighbors x on G n . Thus for every 6" vertex x in G n there is a linear equation 
for — A n u(x); these equations can be stored into a 6" square matrix. The eigenvalues and eigenvec- 
tors of these graph Laplacian matrices are calculated using the eigs command in MATLAB. Table 



6. 1 lists the first twenty eigenvalues for the level-n hexacarpet for n = 7, 



Using our data, we are able to plot the hexacarpet in eigenvalue coordinates. This is analogous 
to harmonic coordinates which are used extensively in the literature on fractals, see |[T8l [201 [33l 
and references therein. Given two eigenf unctions, (p, ip, defined on G n , we plot the ordered pair 
(ip(x),ip(x)) for each x E G n . Figure 6.1 shows the plots of (cpi, (fj) for 2 < i < j < 6 for level 



n = 7. The first eigenfunction <pi is a constant function at each level and is excluded. Figure 6.2 
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(A) (y> 2 ,¥>3) 



(B) ((f2,(pi) 



(C) (</>2,</>0 




(D) ((p2,(p6) 



(E) (<P3,(fi) 



(F) (P3,<P5) 






(G) (¥?3,^6) 



(h) (934,955) 



(i) (y>4,<A>) 



Figure 6.1. Two-dimensional eigenfunction coordinates 



x 



X 



)) of the hexacarpet. 



shows three dimensional plots (</?j, Lfj,Lf k ) for some choices of 2, j, These plots demonstrate that 
the geometry of the hexacarpet is significantly different from the geometry of the triangle, on which 
the hexacarpet construction is based. 

In 0, it is believed that there is a renormalization factor r such that 



A'; 



(6.3) 

where A™ represents the jth eigenvalue of the level n graph Laplacian. According to @, this coeffi- 
cient r is the Laplacian scaling factor (it is denoted R in [6J). This scaling factor satisfies r = Np, 
where iV is the factor that describes how X n grows at each level n or, alternatively, N can also be 
thought of as the number of contraction maps for a carpet, and p is the resistance scaling coefficient 
(see for instance EQjJlE) and references therein). In the case of the hexacarpet, iV = 6 and p is 
not known (in the present state of knowledge it is not possible to compute p theoretically for the 



hexacarpet). Since the renormalized ratios of eigenvalues seems to converge, as seen in table 6.2 
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(a) (^2,^3,^4) 



(b) (^2,^3,^5) 



(c) (ip 2 ,ip 3 ,ipe) 





(d) (^2,^3,^7) 



(e) (^2,^4,^5) 



(F) (^2,^4,<^6) 




(G) {<P2,tp5,<P6) 



(h) (tp 3 ,ip5,<p 6 ) 



(i) (^4,^5,^6) 



Figure 6.2. Three-dimensional eigenfunction coordinates of the hexacarpet. 



we venture Conjecture [T7T]part (1) and (2) in the possibility that a Laplacian and its corresponding 
diffusion process exist on the limit object, the hexacarpet. The best estimate for p will come from 
our estimates for the lowest nonzero eigenvalue, thus giving us an estimate of p m 1.3064. The 
hexacarpet has six contraction mappings, and the natural contraction ratio from Section [2] is 1/2. 
With r = 6p, our calculation above suggests that the spectral dimension is 2 log(6) / log(r) ~ 1.74. 
Part (3) of Conjecture 1 1 . 1 | comes from the fact that the hexacarpet has well defined reflection sym- 
metries, which make it plausible to apply the methods of [|4l [3J . Moreover, we expect logarithmic 
corrections because the approximating graphs have diameters which seem to grow on the order of 
n6 n rather than 6 n , as in Proposition |5.2| and Conjecture |5.4[ There is an extensive literature on 



heat kernel estimates and their relation to functional inequalities (see flU [31 [121 [141 CGI |2T] [301 and 
references therein), which provides general framework and background for our calculations and 
conjectures. 
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Figure 6.3. Shape of the eigenvalue counting function N(X), which is the number 
of eigenvalues less than A, of the 6th level graph for the first 500 eigenvalues. Flat 
intervals correspond to the gaps in the spectrum. 
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1.2801 
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1.2461 


1.3051 


1.3063 


1.3048 


1.3065 


1.3065 


12 




1.2461 


1.3019 


1.3075 


1.3068 


1.3066 


1.3065 


13 




1.1969 


1.6014 


1.0590 


1.3068 


1.3066 


1.3065 


14 




1.1969 


1.2972 


1.3063 


1.3078 


1.3066 


1.3065 


15 




1.2026 


1.3059 


1.3020 


1.3060 


1.3066 


1.3065 


16 




1.2026 


1.2993 


1.3074 


1.3071 


1.3067 


1.3065 


17 




1.1640 


1.3655 


1.2349 


1.3064 


1.3066 


1.3065 


18 




1.1755 


1.4128 


1.2009 


1.3069 


1.3067 


1.3065 


19 




1.1761 


1.5252 


1.1171 


1.3073 


1.3068 


1.3066 


20 




1.1761 


1.2988 


1.3114 


1.3077 


1.3068 


1.3065 



Table 6.2. Hexacarpet estimates for resistance coefficient p given by 



6 A" 
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The reasoning for part (4) of the conjecture is illustrated in Figure [63j as gaps in the spectrum 
correspond to flat intervals in the eigenvalue counting function. It is a worthwhile point of inves- 
tigation because of OTlfTTll . Finally, part (5) of the conjecture is ventured because of recent work 
on spectral zeta functions [8J (see [|22l|32l for background). 
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